The effect of partonic wind on charm quark correlations in high-energy nuclear 

collisions 
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In high-energy collisions, massive heavy quarks are produced back-to-back initially and they 
are sensitive to early dynamical conditions. The strong collective partonic wind from the fast 
expanding quark-gluon plasma created in high-energy nuclear collisions modifies the correlation 
pattern significantly. As a result, the angular correlation function for DD pairs is suppressed at the 
angle A<fi = tt. While the hot and dense medium in collisions at RHIC ( ^/snn = 200 GeV) can 
only smear the initial back-to-back DD correlation, a clear and strong near side DD correlation is 
expected at LHC {y/s]^f = 5500 GeV). 

PACS numbers: 25.75.-q, 12.38.Mh 



The ultimate goal of relativistic heavy ion collision ex- 
periments at RHIC and LHC is to detect the formation 
of a possible new state of matter, the so-called quark- 
gluon plasma (QGP), created in the early stage of the 
collisions. An essential difference between elementary 
particle collisions and nuclear collisions is the develop- 
ment of the collective motion in the latter. The partonic 
collectivity is one of the consequencies of the formation 
of QGP. The collective flow of all hadrons, especially the 
multi-strange hadrons (f> and f2, has been experimentally 
measured [H, 0, 0, 0] ■ These results are consistent with 
the predictions based on hydrodynamic model calcula- 
tions @, 0, ID, 0] and demonstrate strong collective 
expansion at RHIC. However, the issue of local thermal- 
ization of the partonic matter remains to be addressed. 

The observables related to heavy quarks are of par- 
ticular interest in the study of thcrmalization. Unlike 
light quarks which are controlled by chiral dynamics and 
are almost massless in the deconfined phase where the 
chiral symmetry is restored, heavy quarks remain mas- 
sive in QGP and can only be pairly-created in the early 
stage of the collisions at RHIC. In addition, during the 
subsequent evolution, the number of heavy quarks is un- 
changed, because the typical achievable temperature of 
the medium is in general lower than the thresholds for 
heavy quark production. However, the heavy quarks can 
participate in collective motion if, and only if, their in- 
teractions at partonic level occur at a sufficiently high 
frequency. As the light quark and gluon densities are 
much higher than that of the heavy quarks, any collec- 
tive motion of heavy quark hadrons will be a useful tool 
for studying the early thermalization of light quarks in 
high-energy nuclear collisions. 

From lattice calculations [lo| . the thermodynamic 
functions can not reach their Stefan-Boltzmann limits 
even at a temperature of four times the critical value, 
indicating that the system is still strongly coupled at 
very high temperature. Recently, the elliptic flow and 
nuclear modification factors of heavy quarks' decay elec- 



trons are measured [ll|, ll2| at RHIC. In order to explain 
the data, a large drag coefficient in the Langevin equa- 
tion, describing the propagation of charm quarks in the 
strongly coupled QGP, is necessary [HI, [HI]. Two- and 
three-body interactions [HI, [HI] of heavy quarks seem to 
become important as well. Furthermore, from the study 
of J/t/j regeneration in QGP [13], we have found that 
thermalized charm quarks can better describe the J /if) 
averaged transverse momentum square 17, HI- All these 
analysis suggest that heavy quarks have been participat- 
ing in frequent rescatterings at the early partonic stage 
in high-energy nuclear collisions. 

In hard scatterings, heavy quark pairs are produced 
back-to-back at leading order. A typical example for the 
processes is the gluon fusion: gg — > cc. In hadron-hadron 
collisions, although the initial kx kick, high order pro- 
cesses and fragmentation smear the angular correlations, 
the back-to-back peak has been observed in the final state 
DD correlation [19] . In high-energy nuclear collisions, on 
the other hand, a hot and dense partonic medium is ex- 
pected at the early stage. The interactions between the 
charm quarks and the medium will modify the angular 
correlations and even lead to a complete absence of this 
back-to-back correlation [2(| in the final D-meson corre- 
lation function. 

As we emphasized above, the QGP formed in the early 
stage of high-energy nuclear collisions is not a static sys- 
tem, but a fast expanding fireball. The heavy quarks, 
which frequently interact with QGP medium, will expe- 
rience the strong collective flow — the partonic wind. The 
effect of the partonic wind on the DD correlation func- 
tion can be estimated from the transformation between 
two reference frames. Suppose the c and c in a pair move 
radially in QGP and the velocity of the one which prop- 
agates towards the center of the fireball is v c in the rest 
frame of the fluid element. When v c is less than the 
radial flow Vf of the fluid element with respect to the 
laboratory frame, both charm quarks will move together 
towards outside and the initial back-to-back correlation 
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is turned to the same side by the partonic wind! In cen- 
tral Au+Au collisions at RHIC, taking Vf ~ 0.6c 2l|, the 
critical charm quark momentum for the turning of the cc 
correlation is in the order of 1 GeV. Taking into account 
both the partonic wind and the scatterings between the 
heavy quarks and the medium, the turning could occur 
at a much higher momentum. As a result, the partonic 
wind will speed up the disappearance of the DD corre- 
lation and may even change the back-to-back correlation 
to the near side when the wind is strong enough. 

In this Letter, we focus on the mid-rapidity DD az- 
imuth correlation in p+p and A+A collisions at RHIC 
and LHC energies. We first estimate the charm hadron 
correlation for p+p collisions with PYTHIA [22]. The 
motion of charm quarks in the medium is treated as a 
random walk by the Langevin equation [23| . The medium 
itself is described by a 2+1 dimensional Bjorken type hy- 
drodynamic calculation 171, |24j with proper initial con- 
ditions at corresponding collision energies. The effect of 
the partonic wind on the final DD correlation is obtained 
by solving the coupled Langevin equation with the hy- 
drodynamic evolutions. 

We start by reviewing the DD correlation in pp col- 
lisions. The Monte Carlo event generator PYTHIA [22| 
reproduces well the experimentally observed correlation 
measured at fixed target energies [3] . FigQ] shows the 
calculated correlation for p+p collisions at RHIC and 
LHC energies. Compared to the calculation [2(| at lead- 
ing order (LO) pQCD which contains only pair creation 
processes, the next-to-leading order (NLO) pQCD calcu- 
lation includes pair creation with radiative corrections, 
flavor excitation and gluon splitting by switching on the 
initial and final parton showers 25j, |26j . Although this 
treatment is not exactly at next-to-leading order, some 
aspects of the multiple-parton emission phenomenon are 
reproduced 25,26)]. Since cc pairs from flavor excitation 
and gluon splitting do not show explicit angular corre- 
lation, the NLO DD correlation becomes weak at RHIC 
energy and even vanishes at LHC energy, compared to 
the corresponding LO result. At RHIC energy, the cor- 
relation of DD pairs with high transverse momentum pt 
is mostly preserved, since they are mainly from the pair 
production processes. However, at LHC energy, even for 
pairs with p t larger than 5 GeV/c, the back-to-back corre- 
lation is strongly suppressed, since gluon splitting and fla- 
vor excitation become dominant at very high energy [26| . 

Since there is no obvious DD correlation in p+p colli- 
sions at LHC energy, the measurement of smearing this 
correlation in heavy ion collisions becomes meaningless 
at LHC. However, the medium temperature at LHC is so 
high that heavy quarks are much easier to get thermal- 
ized with the medium, and then flow with the medium 
and be pushed by the partonic wind to the same direc- 
tion. Therefore, near side, instead of back-to-back, DD 
correlation might be possible at LHC. 

To explore whether the partonic wind generated in cen- 
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FIG. 1: (Color online) The DD correlation as a function of 
relative azimuth angle A<f> in p+p collisions at sfs = 200 GeV 
(left) and y/s = 5500 GeV (right), calculated by PYTHIA (v. 
6.327) 22]. The transverse momentum depedence is shown 
by the pt cuts (GeV) on single D-mesons. MSEL=4 (thin 
lines) and MSEL=1 (thicker lines) indicate, respectively, the 
calculation to leading order and to next-to-leading order. 



tral relativistic nucleus-nucleus collisions at RHIC and 
LHC influences the DD angular correlation, we employ 
a non-relativistic Langevin approach which describes the 
random walk of charm quarks in a QGP and was first 
described in Ref. 12311 . 



dp/dt = -j{T)p + fj, 



(1) 



where if is a Gaussian noise variable, normalized such 
that {r]i{t)Vj(t')} = Ci(T)5ij5(t ~ t') with i, j indexing di- 
rections. Both the drag coefficient 7 and the momentum- 
space diffusion coefficient a depend on the local temper- 
ature T. We use the same parameterization for 7 as in 
Ref. [H, 7(T) = aT 2 , with a = 2-KT 6 (fm/cT 1 MeV 2 
and also neglect the momentum dependence of 7. a can 
be calculated from the fluctuation-dissipation relation in 
equilibrium, (pf) — a/27, with (p?) = 1.33 m c T and 
charm quark mass m c = 1.5 GeV/c 2 ;23j. 

The above drag coefficient is estimated with pQCD 
adopting a large coupling constant a s = 0.6 21\. The 



recent pQCD calculations [13|, |28( show a factor of 2-3 



smaller drag coefficient, but the quasi- hadronic resonant 
states in QGP at moderate temperature would result in a 
much larger drag coefficient [13j . Since exact values of the 
drag coefficient from lattice QCD are not yet available, 
we will choose in the following numerical calculations the 
drag parameter a in a reasonable region. The compari- 



son of a similar Langevin calculations 1J| for Raa and 
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V2 of non-photonic electrons with recent PHENIX data 
shows that the data favor a large drag coefficient 
We checked that the drag parameter a corresponding to 
the diffusion coefficient D = 3/(2itT) in [H], is about 
6 • 1(T 6 (fm/c)- 1 MeV~ 2 . 

For simplicity, we focus on charm quarks generated 
at mid-rapidity in central collisions. We assume that 
the produced partonic plasma reaches local equilibrium 
at time r n . After that, the plasma evolves according to 
the 2+1 dimensional Bjorken's hydrodynamics which de- 
scribes the collective flow of the bulk matter very well at 
RHIC Hi, 



d T E + V ■ M 

d T M x + V ■ (M s v) 
V 



d T M, 



-(E + p)/t , 

-m x /t - d xP 

-M v /t— d y p 
d T R + V • (Rv) = -R/t 



(M y v) 



(2) 



with the definitions E — (e + p)j — p, M = (e + p)j v 
and R = 777,, where 7 is the Lorentz factor, and e,p and 
v are the energy density, pressure and transverse velocity 
of QGP. To close the hydrodynamic equations, we take 
the equation of state of ideal quark-gluon plasma [3ll ] 
with first order phase transition at T c = 165 MeV. The 
initial condition for central collisions of heavy nuclei at 
RHIC is the same as in 24[. For Pb+Pb collisions at 



LHC energy, the initial baryon density is neglected and 
the maximum initial entropy density is assumed to be 5.5 
times the value (about 110 fm -3 [32j) at RHIC. Taking 
into account the earlier equilibrium time at LHC, this 
initial value is reasonable according to the multiplicity 
prediction at LHC 33]. For a simple comparison, the 
maximum temperature To is 340 MeV at To = 0.6 fm/c 
at RHIC and 610 MeV at r = 0.3 fm/c at LHC. 

Before the time To, the plasma changes rapidly and is 
not fully equilibrated. To simplify the calculation, wc 
assume, as done in [23| , that the plasma is already in 
equilibrium before tq with a constant temperature T = 
To, and the charmed quarks diffuse in the same way as 
they do after to- 

The initial charm quarks are generated with PYTHIA 
at a time of the order of l/m c ~ 0.1 fm/c. In order to 
isolate the effects purely due to parton-parton rescatter- 
ing in the outlined medium, we use a delta function for 
fragmenting a charm quark into a charmed hadron at the 
hadronization stage. 

The radius r where each pair is created is randomly 
generated from a distribution reflecting the number of 
binary nucleon-nucleon collisions taking place at that ra- 
dius, p(r) oc (R A — r 2 )2nr, where Ra is the geometric 
radius of the colliding nuclei. 

The evolution of the charm momentum with the 
Langevin equation ends when T reaches the critical tem- 
perature T c or when the charm quark leaves the QGP 
volume. To get a first estimation of the QGP effect on 



the pair correlation, we omit the contribution from the 
mixed phase and hadronic phase. 

It is necessary to note that the hydrodynamic equa- 
tions determine not only the local temperature T which 
characterizes the drag coefficient but also the local fluid 
velocity v which controls the partonic wind effect on the 
cc correlation. Since the Langevin equation is defined 
in the local rest frame of fluid element, we should make 
a transformation from the rest frame to the laboratory 
frame which is governed by the fluid velocity v. 
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FIG. 2: (Color online) The DD correlation as a function of 
relative azimuth angle A(j> with different drag parameter a 
for central Au+Au collisions at RHIC (left) and Pb+Pb col- 
lisions at LHC (right). The unit of the drag coefficient a is 
[10~ 6 (fm/c)- 1 MeV" 2 ]. In the PYTHIA calculation to next- 
to-leading order, no pt cut was applied. The initial tempera- 
ture To and thermalization time to are, respectively, 340 MeV 
and 0.6 fm/c at RHIC and 610 MeV and 0.3 fm/c at LHC. 

Fig. [2] shows the results of DD angular correlation for 
central y/s = 200A GeV Au+Au collisions at RHIC (left 
plot) and V« = 5500A GeV Pb+Pb collisions at LHC 
(right plot). The thin solid lines with a — are the 
corresponding results directly from PYTHIA with par- 
ton showers turned on. At RHIC, with increasing drag 
coefficient the initial back-to-back DD correlation grad- 
ually become more and more flat. At a ~ 6 ■ 10 -6 
(fm/c) _1 MeV -2 , the correlation is totally washed out, 
see left plot. Note that the recent PHENIX data [HI 
seem to favor the relatively large value of the drag coef- 
ficient. The experimental measurement of this smearing 
effect in the heavy quark correlation, especially for those 
pairs at high p tl can shed a light on the thermalization 
status of the medium at RHIC [2(J. At LHC, one starts 
with an almost flat correlation, see the thin solid line in 
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the right plot. The near side DD correlation becomes 
more and more visible with increasing the drag coeffi- 
cient. For a ~ 4 ■ 10~ 6 (fm/c) _1 MeV~ 2 , the near side DD 
correlation is already very strong. The larger the drag co- 
efficient means the stronger the interactions. This result 
indicates that cc pairs are thermalizcd with the medium 
quickly and they are moving in the same direction as the 
partonic wind. 

Note that there are more than one pair of DD pro- 
duced in high energy nuclear collisions. At mid-rapidity 
\y\ < 1) the total number of charm pairs is in the order of 
2-4 and 6-10 at RHIC and LHC, respectively. The mul- 
tiple production of un-correlated D and D mesons will 
give an additional background to the measured angular 
correlations. Our tests show that this background can be 
removed by the mixed event method. The mixed-event 
subtracted correlation function should be compared with 
the calculated results shown in Fig. 2. 

In summary, we investigated the partonic wind effect 
on DD angular correlation in high-energy nuclear colli- 
sions. At RHIC, while there is still sizable DD back-to- 
back correlation in p+p collisions, it is washed out in the 
most central Au+Au collisions. At LHC, the pQCD high 
order contributions to c-quark production become dom- 
inant and there is no more clear DD angular correlation 
even in elementary collisions. However, charm quarks 
are thermalized with the medium quickly at LHC. As a 
result, a near side DD correlation, driven by the strong 
partonic wind, will replace the flat angular distributions. 
Such dramatic changes in the DD angular correlations, 
originated from the strong interaction between charm 
quarks and the hot medium, reflected by the large drag 
coefficient (> 4 • 10~ 6 (fm/c) _1 MeV~ 2 ) in the Langevin 
equation in our treatment, should be considered as a sig- 
nature of the strongly coupled QGP formation at RHIC 
and LHC. 
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